Magnetic Ordering in Carbon Nanotubes 
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The physical origin of the mechanism of the formation of ferromagnetic ordering in carbon nan- 
otubes (NTs) produced by high energy ion beam modification of diamond single crystals in (110) 
and (111) directions has been found. It is determined by asymmetry of spin density distribution 
of Su-Schrieffer-Heeger topological soliton lattice formed in ID Fermi quantum liquid state of n- 
electronic subsystem of given NTs. The phenomenon of formation of antiferromagnetic ordering 
coexisting with superconductivity at room temperature in carbon NTs, produced by high energy 
ion beam modification of diamond single crystals in (100) direction is discussed. 
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I. INTRODUCTION AND BACKGROUND 

It is well known that all substances on the whole are 
magnetics. At the same time, it is also well known, that 
classical magnetic ordering is existing in the substances, 
which arc built from the atoms with unfilled inner atomic 
d- or /-shells or include given atoms in their elementary 
units. In other words, it is the substances, elementary 
units of which include transition chemical elements with 
unfilled atomic 3d-, 4d-, 5d-, 6d-shells, or 4f, 5f-shells 
of rare earth elements. Carbon does not refer to given 
group. Nevertheless, there are at present a number of 
communications on magnetic ordering in carbon and car- 
bon based materials. 

On the experimental revealing of magnetic ordering 
in carbon structurally ordered systems was reported for 
the first time during the IBMM-Conference in Knoxville, 
TN, USA [1] and it was confirmed in report on E-MRS 
Conference in Strasbourg, France 0. Let us remark, 
that the first report almost in the same time on mag- 
netic ordering in structurally non-ordered carbon mate- 
rials is the work [3[ , where ferromagnetic ordering in py- 
rolytic carbon, produced by chemical vapour deposition 
(CVD) method using adamantane to be source material, 
was found. Let us also remark, that simultaneously, the 
reports 0, were the first reports on the formation 
by high energy ion beam modification (HEIBM) of di- 
amond single crystals structurally and magnetically or- 
dered quasi-one-dimensional (quasi-lD) system along ion 
tracks, that is new carbon allotropic form, which was 
identified with nanotubes NTs, incorporated in diamond 
matrix in direction, precisely coinciding with ion beam 
direction. Given system possesses by a number of inter- 
esting physical properties, reported in 0, 0, 0- When 
concerne the magnetic ordering, it was established from 
the study of temperature dependence of electron spin res- 
onance (ESR) absorption intensity, that, for instance, in- 
corporated nanotubes, produced by neon HEIBM of di- 
amond single crystal along (100) crystallographic direc- 
tion, possess by weak antiferromagnetic ordering 0], 0, 
0. At the same time, copper HEIBM with implantation 



direction along (111) crystal axis, nickel HEIBM with 
implantation direction along (110) axis 0, 0, an d 
boron HEIBM of polycrystalline diamond films with im- 
plantation direction transversely to film surface [5| lead 
to formation of NTs, incorporated in diamond matrix, 
which are possessed by ferromagnetic ordering. It was 
established directly by observation of ferromagnetic spin 
wave resonance (FMSWR) 0, 0, 0. The first obser- 
vations of FMSWR in carbon materials and in the ma- 
terials, which are not traditional magnetic materials, on 
the whole, were reported in above cited works. It was 
established, that magnetic ordering is inherent property 
for carbon electronic system and it is not connected with 
magnetic impurities, since starting samples were selected 
in that way, that the absolute spin number of paramag- 
netic impurities and paramagnetic structural imperfec- 
tions at all did not exceed the value 10 12 spins. Very 
recently 0, antiferroelectric ordering has been found in 
the same pure carbon allotropic form - quasi-lD carbon 
zigzag-shaped nanotubes (CZSNTs), obtained by boron- 
and copper-HEIBM of diamond single crystals in (In- 
direction. It was established by means of the detection 
of new optical phenomenon - antiferroelectric spin wave 
resonance ( AFESWR) , which was theoretically described 
and experinmentally confirmed for the first time by in- 
frared (IR) spectroscopy study of carbynes in It 
seems to be very new property of pure carbon allotropic 
forms - quasi-lD CZSNTs and carbynes. Moreover, on 
the observation of antiferroelectric ordering has been re- 
ported in [ll| for the first time (to our knowlegde) for 
all matter systems on the whole. Given results mean, 
that pure carbon in the form of quasi-lD CZSNTs and 
carbynes is multiferroic system. 

Let us remark, that experimental observation of mul- 
tiferroicity in quasi-lD CZSNTs and carbynes means the 
breakdown of space inversion symmetry along CZSNT 
hypercomplex symmetry axis and carbyne chain symme- 
try axis. In the case of CZSNTs, it agrees well with 
the model of quasi-lD CZSNTs Q, @, based on bond 
dimerization in all chain components of quasi-lD CZSNT 
along its hypercomplex symmetry axis z, which actually 
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leads to inversion symmetry breakdown along given axis. 
Therefore, the experimental observation of antiferroele- 
cricity of quasi-lD CZSNTs, necessary condition of which 
is the evident prediction of the model (appearance of 
nonzero polarisation by atomic displacements) , proposed 
in [H, can be considered to be additional argument in 
favour of given model. 

The aim of given work is to study in more details 
the properties of non-cylindrical nanotubes possessing in- 
stead Coo symmetry axis the only C4 symmetry axis and 
to establish the mechanisms of formation of magnetic and 
electric ordering in NTs both with Coo symmetry axis 
and with C4 symmetry axis, produced in diamond single 
crystals by high energy ion implantation. 



II. EXPERIMENTAL TECHNIQUE AND 
RESULTS 



Samples of type Ha natural diamond, implanted by 
high energy ions of nickel (the energy of ions in ion beam 
was 335 MeV) have been studied. Paramagnetically pure 
samples have been selected so that the absolute spin num- 
ber did not exceed the value rj 10 12 spins in each of the 
sample used before implantation. Ion implantation was 
performed along (100) crystal direction (ion beam dose 
was 5 x 10 13 cm -2 ), and in (110) crystal direction (ion 
beam dose was 5 x 10 14 cm -2 ) transversely to sample 
(100) and (110) planes correspondingly uniformly along 
all the plane surface. The temperature of the samples 
during the implantation was controlled and it did not 
exceed 400 K. ESR spectra were registered on X-band 
ESR-spectrometer "Radiopan" at room temperature by 
using of TE102 mode rectangular cavity. The ruby stan- 
dard sample was permanently placed in the cavity on its 
sidewall. One of the lines of ESR absorption by Cr 3+ 
point paramagnetic centers (PC) was used for the cor- 
rect relative intensity measurements of ESR absorption, 
for the calibration of the amplitude value of magnetic 
component of the microwave field and for precise phase 
tuning of microwave field. It was possible owing to unsat- 
urating behavior of ESR absorption in ruby in the range 
of the microwave power applied, which was w 100 mW in 
the absence of attenuation. Unsaturable character of the 
absorption in a ruby standard was confirmed by means 
of the measurements of the absorption intensities in two 
identical ruby samples in dependence on the microwave 
power level. The first sample was standard sample, per- 
manently placed in the cavity, the second sample was 
placed in the cavity so that the resonance line intensity 
was about 0.1 of the intensity of corresponding line of 
the first sample. They were registered simultaneously 
but their absorption lines were not overlapped owing to 
slightly different sample orientations. The foregoing in- 
tensity ratio was preserved for all microwave power val- 
ues in the range used, which indicates, that really ruby 
samples are good standard samples in ESR spectroscopy 
studies. 
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Figure 1: Angular dependence of effective g- value for the third 
FMSWR-mode, observed in NTs, incorporated in diamond 
single crystal by nickel ion beam direction transversely (110) 
sample plane 



The studies of FMSWR in the sample implanted along 
(110) crystal direction were described in Q, however the 
mechanism of the ferromagnetically ordered state was not 
proposed. Let us summarise the results, which seems to 
be essential to establish the detailed processes, leading to 
ferromagnetism of incorporated carbon NTs, produced 
by (110) implantation. The strong anisotropy of spin 
wave excitation have been observed in unannealed sam- 
ples. FMSWR was registered the only in the case, when 
the implantation plane was perpendicular to the vector 
Hi of magnetic component of the microwave field. The 
intensity of FMSWR absorption was increasing with tem- 
perature of isochronal 20-minutes annealing up to ~ 700 
Celsius degrees. It is interesting, that after the anneal- 
ing at 350 Celsius degrees and by higher temperatures 
FMSWR absorption was detected also by the vector Hi 
of magnetic component of the microwave field position, 
being to be parallel to sample implantation plane. The 
parameters of FMSWR absorption, that is splitting pa- 
rameter, which can be characterised by, for instance, ef- 
fective g-values of its modes, in particular by g-value of 
the third mode, Figure 1, and the modes' linewidths and 
intensities were strongly anisotropic. Angular dependen- 
cies of the linewidth, FMSWR absorption amplitude and 
intensity of the third FMSWR-mode are presented in Fig- 
ures 2 to 4. 

The splitting parameter anisotropy seems to be non- 
trivial, it characterises by two maxima - at 54.7 degrees, 
that is, in [111] direction of diamond lattice, and in di- 
rection about 17 degrees from [100] crystal direction. It 
is interesting, that unlike to (111) ion beam modifica- 
tion [|| , the angular dependence of linewidth has the 
resemblance with the angular dependence of g-value, it 
characterises also by two maxima at the same angles, 
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Figure 2: Angular dependence of the linewidth of the third 
FMSWR-mode, observed in NTs, incorporated in diamond 
single crystal by nickel ion beam direction transversely (110) 
sample plane 
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Figure 3: Angular dependence of FMSWR absorption ampli- 
tude of the the third FMSWR-mode, observed in NTs, incor- 
porated in diamond single crystal by nickel ion beam direction 
transversely (110) sample plane 



see Figure 2. (At (111) ion beam modification minimum 
of g-value, for instance, for main FMSWR-mode corre- 
sponds to maximum of linewidth value) . Minimal effec- 
tive g-value and minimal value of linewidth for the third 
mode are achieved at [110] direction and they are equal 
to 2.1526 and to 100 G respectively, maximal effective 
g-value and maximal value of linewidth correspond to 
[111] direction of diamond lattice and are equal respec- 
tively to 2.3845 and 360 G. Coinciding of the directions 
of anisotropy maxima with diamond lattice characteris- 



Figure 4: Angular dependence of the FMSWR absorption 
intensity of the third FMSWR-mode, observed in NTs, incor- 
porated in diamond single crystal by nickel ion beam direction 
transversely (110) sample plane 
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Figure 5: Spectral distribution of ESR absorption intensity 
in diamond single crystal, implanted by high energy nickel 
ions by beam direction transversely (100) sample plane, the 
sample was rotated in (011) plane, Ho\\ [100] crystal axis 



tic directions is indication of the role of host diamond 
lattice on the characteristics of magnetic ordering. It 
is also strong evidence, that magnetic ordering is de- 
termined by electron-electron correlations in valence 7r- 
subsystem of electronic system (but not in er-subsystem), 
since the only 7r-subsystem can be subjected to so strong 
response on the surrounding lattice presence by NT for- 
mation, which is starting at more lower ion fluences, and 
strong response seems to be preserved by entire HEIBM 
of near surface region in diamond single crystal. Ampli- 
tude dependence, Figure3, is characterised by one clearly 
pronounced maximum near 35.3 degrees direction of di- 
amond lattice, that is, it coincides also with symmetry 
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Figure 6: Spectral distribution of ESR absorption intensity 
in diamond single crystal, implanted by high energy nickel 
ions by beam direction transversely (100) sample plane, the 
sample was rotated in (011) plane, _Ho|| [HI] crystal axis 



Figure 8: Dependence of absorption amplitude of the left line 
in ESR spectrum of NTs incorporated in diamond single crys- 
tal on magnetic component of microwave field at Ho\\ [100] 
crystal axis 




Figure 7: Spectral distribution of ESR absorption intensity 
in diamond single crystal, implanted by high energy nickel 
ions by beam direction transversely (100) sample plane, the 
sample was rotated in (011) plane, the angle between Ho and 
[100] crystal axis is 60 degrees 



Figure 9: Dependence of absorption amplitude of the right 
broad line in ESR spectrum of NTs incorporated in diamond 
single crystal on magnetic component of microwave field at 
Ho 1 1 [100] crystal axis 



direction of diamond lattice. It is seen, that amplitude 
dependence is also strongly anisotropic. So, amplitude 
change is equal to 4.8 times. It is additional indication 
that mechanism of magnetic ordering is connected with 
the formation of ordered structures just in 7r-subsystem 
of incorporated NTs. Simultaneously, it is the display 
of the fact, that anisotropy of FMSWR absorption by 
magnetically ordered NT-structure and the anisotropy 
of magnetically ordered NT-structure itself are different. 



Two not very pronounced maxima at 54.7 degrees and 
in direction near 17 degrees from [100] crystal are ap- 
peared additionally to the maximum near 35.3 degrees 
direction in the angular dependence of the intensity, Fig- 
ure 4. They are determined by contribution of linewidth 
anisotropy, compare Figure 4 and Figure 2. 

The ESR spectra observed in carbon nanotubes (NTs) , 
produced by nickel high energy (100) ion beam modifi- 
cation of natural diamond single crystals, are presented 



5 



in Figures 5 to 7 in crystal directions [100], [111] and 60 
degrees from [100] correspondingly. The line in the range 
(1865 - 1981) G (it is indicated by arrows in Figure 4) 
is the absorption line by ruby standard, it isshifted to 
botton in Figure 5 and it is removed in the same range 
in Figure 6. The most intensive lines, belonging to the 
sample studied, were excited spontaneously the only by 
very precise orientation of external static magnetic field 
Ho in the direction coinciding with implantation direc- 
tion and resulting spectrum was consisting of three lines, 
at that two lines have rather large anisotropic linewidths. 
Let us designate given lines by Rb for the right broad line 
and by L for the left line. The right broad line was over- 
lapped with relatively narrow almost isotropic line, des- 
ignated by R n (given line was observed by usual sample 
orientation). Additionally, very broad strongly intensive 
anisotropic absorption with two dip positions at ~ 2410 
G and ~ 2892 G in the direction corresponding [111] di- 
amond lattice direction, Figure 6, is detected, which are 
overlapped by static magnetic field direction in 60 de- 
grees from [100] diamond crystal direction, Figure 7. It 
seems to be the display of the fact, that the symmetry of 
the interaction, leading to the appearance of very strong 
absorption is determined by inherent symmetry of NTs, 
produced by [100] HEIBM, which is not connected with 
potential effect of diamond lattice presence. 

Dependence of absorption amplitude of L-line and Rb 
line on magnetic component of microwave field at fixed 
orientation of polarising magnetic field i?o|| [100] crystal 
axis have been studied, Figures 8 and 9. It is seen from 
comparison of the Figures 8 and 9, that given dependen- 
cies are quite different. The dependence, presented for Li- 
line in Figure 8 is superlinear. It is similar to the depen- 
dencies, which were earlier observed in the samples, mod- 
ified by HEIBM with copper, neon, nickel ions (however 
with dose 5 x 10 14 ) [4| , [7J , [6| , that is, in the case of entire 
near surface layer modification. In the studied sample 
(integral dose is 5 x 10 13 ), individual NTs are isolated by 
diamond structure, nevertheless the superlinear depen- 
dence is taking place, which seems to be unexpected. Let 
us remark, that the initial part of the curve, presented in 
Figure 8, can in the principle be approximated by a linear 
dependence (dashed line) , although it is clear from com- 
parison with the approximation of the whole curve, solid 
line in Figure 8, that even initial part, strongly speaking, 
is not linear. Solid line in Figure 8 is the polinomial fit 
with the function f(x) = 6rj + b\x + b 2 x 2 + 63a; 3 + 64X , 
where b = -0.17117, 61 = 208.92305, b 2 = -139.06624, 
63 = 159.14424, 64 = -33.90983. 

Dependence of absorption amplitude of the right broad 
line Rb in ESR spectrum of NTs on magnetic component 
of microwave field is strongly nonlinear. It is charac- 
terised for the values of relative magnetic component of 
microwave field Hi/H^ in the range (0-0.75) by usual 
saturating law, but in the range (0.75-1) it acquires 
prominent superlinear nonsaturating character. The de- 
pendence for ESR absorption kinetics in the form of pre- 
sented in Figure 9 is observed in ESR-spectroscopy for 
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Figure 10: Angular dependence of g-factor of the left line in 
ESR spectrum of NTs incorporated in diamond single crystal, 
the sample was rotated in (011) plane 
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Figure 11: Angular dependence of ESR absorption intensity of 
the left line in ESR spectrum of NTs incorporated in diamond 
single crystal, the sample was rotated in (011) plane 



the first time. It can be approximated by the solid line, 
which represents itself the polynomial fit in accordance 
with relation Y(x) = A + B x x + B 2 x 2 + B 3 x 3 + B 4 x 4 , 
where A = 1.14618, B x = 0.77956, B 2 = -0.00159, 
B 3 = -3.03868£ - 5, B A = 1.40072.B - 7. Angular 
dependence of g-factor of the left line L in ESR spec- 
trum of NTs in the sample studied is presented in Figure 
10. It consist of two branches. One branch is in the an- 
gle range 0-60 degrees from [100] crystal lattice direction 
(which is coinciding with NT axis direction), the second 
branch is in the angle range 60-90 degrees. Let us re- 
mark that the connection point of two branches, equaled 
to 60 degrees for the g- values of L-line is coinciding with 
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Figure 12: Angular dependence of linewidth of the left line in 
ESR spectrum of NTs incorporated in diamond single crystal, 
the sample was rotated in (011) plane 



Figure 14: Angular dependence of the cavity Q-factor with 
the sample implanted by high energy nickel ions by beam 
direction transversely (100) sample plane, the sample was ro- 
tated in (011) plane 
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Figure 13: Angular dependence of asymmetry extent A/B of 
the left line in ESR spectrum of NTs incorporated in diamond 
single crystal, the sample was rotated in (011) plane 



the point of the junction of two dips in the very broad 
(and consequently very intensive) absorption, testifying 
on the same or related nature of the resonance processes, 
which are responsible for the appearance of L-line and 
very broad lines. The deviation of g- values from free 
electron value g = 2.0023 is very large, at that the min- 
imal value is achieved in the range 16-20 degrees from 
the [011] direction in diamond lattice and it is equal to 
sa 2.0719, maximal g-value corresponds to NT axis di- 
rection, that is to [100] crystal lattice direction and it 
is equal to « 2.3120. Given values are characteristic for 



the systems with the strong magnetic ordering. Conse- 
quently, we have obtained the direct proof of the sponta- 
neous transition of NT system, incorporated in diamond 
lattice, in the state with the strong magnetic ordering. 
Angular dependence of ESR absorption intensity of the 
left line L has qualitatively opposite character. The max- 
imal absorption value corresponds to the direction, being 
to be transversal to NT axis, which is coincides with [011] 
direction in diamond lattice, Figure 11. Additional max- 
imum is observed at 60 degrees from given direction. Let 
us remark that both the maxima in angular dependence 
of ESR absorption intensity of the line L are observed 
also in angular dependence of linewidth, Figure 12, in- 
dicating that the main features in given dependence are 
governed by angular dependence of linewidth. It is con- 
firmed also by that, the very pronounced maximum in 
the angular dependence of absorption amplitude of line 
L is in the range 10-20 degrees from [100] crystal lattice 
direction (the corresponding Figure is not presented). At 
the same time it is not pronounced in angular dependence 
of absorption intensity, see Figure 11. Especially inter- 
esting, that the line L is asymmetric, Figure 13. However 
numerical values of the ratio A/B of the asymmetry ex- 
tent are quite different in comparison with those ones by 
usual Dyson shape [24]]. The value A/B is equal 2.55 for 
static (immobile) PC in conductive media and it is deter- 
mined by by the space dispersion contribution p25j , which 
is appeared in conductive media. It corresponds to the 
space dispersion contribution and absorption contribu- 
tion to resulting ESR response equal to (1 : 1) [25[. The 
value A/B for absorption derivative is increasing from 
2.55 to more than 19 for mobile PC in dependence on 
velocity of the spin diffusion (26|. It is seen from Fig- 
ure 13 that the ratio A/B is less than 1 and it is strongly 
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anisotropic. The maximal A/B value is near [111] crystal 
lattice direction and it is equal to 0.83, Figure 13. The 
minimal A/B value is near 30 degree from [100] crystal 
lattice direction and it is equal to 0.49, Figure 13. The 
presence of the ruby standard has allowed to control the 
cavity Q-factor, Figure 14. It seems to be substantial 
that Q-factor is increasing in the range where A/B value 
is decreasing, that is near 30 degrees and in the range 
near 70-80 degrees, compare Figure 14 and Figure 13. 



III. DISCUSSION 

Although it was established, that copper HEIBM with 
implantation direction along (111) crystal axis, nickel 
HEIBM with implantation direction along (110) axis, 
and boron HEIBM of polycrystalline diamond films with 
implantation direction transversely to film surface lead 
to formation of NTs, incorporated in diamond matrix, 
which possess by ferromagnetic orderingn Q, Q, 
Q the nature of given state and mechanism, leading to 
its formation was not theoretically described. To solve 
given task, it seems to be necessary to know the nature 
of charge and spin carriers and the mechanisms of car- 
rier transport and interactions of charge and spin carri- 
ers both between themselves and with phonons and pho- 
tons. There seems to be paramount significant the same 
task for nanoelectronics, spintronics and for the other 
branches of nanotechnology. There is existing in the the- 
ory of ID electronic systems in particular in the theory of 
conducting NTs the following concept, which was start- 
ing with the work of Tomonaga in 1950 [l6| and with the 
work by Luttinger in 1963 [17[, when it has become clear 
that the electron-electron interaction destroys the sharp 
Fermi surface and leads to a breakdown of the Landau 
Fermi liquid (LFL) theory. The resulting non-LFL state 
is commonly called Luttinger liquid (LL), or sometimes 
Tomonaga-Luttingcr liquid to describe the universal low- 
energy properties of ID conductors. LL behavior is char- 
acterized by pronounced power-law suppression of the 
transport current and the density of states, and by effect 
of spin-charge separation. The nature of the spin and 
charge carriers was also established. They are chargeless 
spin 1/2 quasiparticles - spinons and spinless quasipar- 
ticles with the charge ± e - holons. The universality of 
LL description means that the physical properties do not 
depend on details of the model, the interaction poten- 
tial, and so on, but instead are only characterized by a 
few parameters - critical exponents. Quite remarkably, 
the LL concept is believed to hold for arbitrary statisti- 
cal properties of the particles, that is both for fermions 
and bosons. It provided a paradigm for non- Fermi liquid 
physics. It is argued in many works, that the single- 
wall carbon nanotubes (SWCNTs), considered to be ID 
objects (it is not always correct, especially for standard 
NTs with diameter in several nanometers) can be de- 
scribed the only in the frame of LL concept. Moreover, 
SWCNTs are considered to be the best model system of 



the LL state demonstration. Really, power-law behavior 
was observed experimentally by measuring the tunneling 
conductance of SWNTs in dependence on temperature 
and voltage. Electron force microscopic measurements 
showed also the ballistic nature of transport in conduct- 
ing SWNTs. At the same time spin-charge separation 
has not been observed so far. Therefore, the existing 
viewpoint seems to be insufficiently grounded. 

It has to be remarked, that both the models LL and 
LFL are the models of ideal (and even oversimplified) 
quantum liquids, since they do not take into account the 
nonlinearity of the fcrmion spectrum on the one hand and 
electron-phonon interactions on the other hand. In fact 
both the models describe not strongly adequately the real 
processes, since the changes in charge state of arbitrary 
atom in ID chain to be the result of electron-electron 
interaction are always accompanied by the changes in 
phonon subsystem (and vice versa). It is consequence of 
generic coupling between operators of creation and anni- 
hilation in electron subsystem and in phonon field. Let 
us also remember, that key argument for insertion of the 
notion " Luttinger liquid" is in fact the simplification, de- 
termined by linearization of the generic spectrum of par- 
ticles in neighborhood of Fermi points in k-space. Just 
given simplification has led to divergencies arising in the 
perturbation theory in ID-case and it does not means 
that ID Fermi liquid description is incorrect in general 
case. Consequently, the description of NTs the only in 
the frames of LL concept seems to be also oversimplifi- 
cation. Moreover, it is showed in [l(| that the concept 
of description of ID correlated electronic systems in the 
frame of ID Fermi liquid (FL) can be renewed. The most 
interesting it was found, that the Fermi liquid concept 
can be applied just to ID carbon NTs, while 2D NTs 
and graphene can be in principle described in the low en- 
ergy physics processes in the frame of LL model. It was 
considered in [l(| the concept of ID FL on the example of 
well known ID system - irans-polyacetylene (t-PL), that 
is, it is in fact the gene ralization of well known Fermi 
gas SSH- model [HI, [HI, which, however, in distinction 
from LFL and LL models, takes into account the electron- 
phonon interaction. The subsequent generalization, for 
instance, for application of given model immediately to 
quasi-lD carbon zigzag shaped nanotubes can be easily 
obtained by using of hypercomplex number theory like 
to the works @, where hypercomplex number the- 
ory was applied for the interpretation of quantum optics 
effects in NTs. 

Let us summarise briefly for the convenience of the 
readers the results, presented in (lfj. 

The Born-Oppcnheimer approximation was considered 
and starting Hamiltonian is the following 

H(u) = H {u) + H v ,t( u ) + (1) 
The first term in ((TJ) is 
P 2 

Ho{v) = ^2^2(i^a+ t3 a m ,a + Ku 2 m a+ iS a m ,a)- (2) 
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It represents itself the sum of operator of kinetic energy 
of CH-group motion (the first term in $2$) and the op- 
erator of the cr-bonding energy (the second term). Co- 
efficient K in @ is effective er-bonds spring constant, 
M* is total mass of CH-group, u m is configuration coor- 
dinate for m-th CH-group, which corresponds to trans- 
lation of m-th CH-group along the symmetry axis z of 
the chain, m = 1,N, N is number of CH-groups in the 
chain, P m is operator of impulse, conjugated to configu- 
ration coordinate u m , m = 1, N, s , a m s are creation 
and annihilation operators of creation or annihilation of 
quasiparticle with spin projection s on the m-th chain 
site in cr-subsystem of t-PA. The second term in ([1]) can 
be represented in the form of two components and it is 



#7r,t(«) = #7r,i ( u ) +Ti7T.t, ai (u) = 

EE^^+M 5 ™^ + C+ jS C m+ l jS )] + (3) 

771 S 

(-l) m 2a lU )(c+ +ls c„ hs + c+ s c m+hs )}, 



where c+ s , c myS are creation and annihilation operators 
of creation or annihilation of quasiparticle with spin pro- 
jection s on the m-th chain site in 7r-subsystem of t-PA. 
It is the resonance interaction (hopping interaction in 
tight-binding model approximation) of quasiparticles in 
7r-subsystem of t-PA electronic system, which is consid- 
ered to be Fermi liquid, and in which the only constant 
and linear terms in Taylor series expansion of resonance 
integral about the dimerized state are taking into ac- 
count. 

Operator "H(u) is invariant under spatial translations 
with period 2a, where a is projection of spacing between 
two adjacent CH-groups in undimerized lattice on chain 
axis direction, and which is equal to 1.22 A. It means, 
that all various wave vectors k in fc-space will be in re- 
duced zone with module of k in the range —■£-< k < 

t 1 ° 2a — — 2a 

[151 ] - It can be considered like to usual semiconductors 
to be consisting of two subzones - conduction (c) band 
and valence (v) band. Then it seems to be convenient to 
represent the operators {c+ s }, {c m . s }, m = 1, N, in the 
form 



{cm,s} — {c m s } + {c m s }, 



(4) 



related to tt — c- and tt — u-band correspondingly, and to 
define /c-space operators 



{eg} = {-=^^(-l) m+1 exp(- i fcm a )cWj, 

* m s 

{eg} = {^£E exp(-ifema)cW }, 



(5) 



m = l,N. 

The cr-operators {a+ s } and {a ms }, m = 1, N can also 
be represented in the form like to Q for 7r-operators and 



analogous to ([5]), transforms can be defined. Then the 
expression for the operator T-Lq{u) can be rewritten 



P 2 



2M* 



(7,C-(7,C . ^ -\-a,v (J,V\ 

a k \ +a h ' a,' 



Ku 2 m )x 



(6) 



where a k s , a k s and a k s , a k s are cr-operators of cre- 
ation and annihilation, related to tr-c-band and to a-v- 
band correspondingly. The independence of \u m \ on m, 

P 2 

m = 1,N, means, that the expression y^jp 

independent on m. Then we obtain 



Kut 



«o(«) = E Efe + + O 

k s 



(7) 



where nZ' c s and nZ' v are operators of number of a- 
quasiparticles in tr-c-band and c-w-band correspondingly. 
The expression for i-L-„,t a {u) in terms of {eg} and 

{e g} is coinciding with known corresponding expression 
in 111, QJ and it is 



Hnfy ( u ) = E E 2t ° cos ka ^l 



+(c)~(c) 
c k,s 



'fc,s L fe,si (8) 



The expression for the second part of operator ?^ T) t(?x) 



in terms of {eg} and {ci^} is also coinciding in its form 
with known corresponding expression in [14j . [TEj and it 
is given by 



s(«) 



'ft-K^ai ( u ) = E E 4 a i u sin feaC 



i+(«)^(c) 



*+(c).(o)n 



(9) 



where subscript ai in Hamiltonian designation indicates 
on the taking into account the part of electron-phonon 
interaction, connected with resonance interaction (hop- 
ping) processes. The expression for the operator i-L^ tV ,[u), 
which describes the part of electron-phonon interac- 
tion, determined by interaction between quasiparticles 
in Fermi liquid state of 7r-subsystem in terms of {eg} 
and {eg} is the following 

n^,u, a2 {u) = E E E "2^' fc/ ' s )4!.4!? g M g M- 

k k' s 

(10) 

Here the consideration is restricted by the taking into 
account the contribution of the term, corresponding to 
interaction between the quasiparticles in different bands, 
which seems to be the most essential. Physically the iden- 
tification of linear on displacement u part of both reso- 
nance interaction (hopping) and the pairwise interaction 
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of quasiparticles in 7r-subsystem between themselves with 
electron-phonon interaction is understandable, if to take 
into account, that by atomic CH group displacements the 
phonons are generated, which in its turn can by release 
of the place on, for instance, (CH) m group, to deliver 
the energy and impulse, which are necessary for transfer 
of the quasiparticle (electron) from adjacent (to — 1)- or 
(to+ Imposition in chain in the case of resonance interac- 
tion (hopping). For the case the pairwise interaction of 
quasiparticles, it means, that its linear on displacement u 
part is realized by means of phonon field, which transfers 
the energy and impulse from one quasiparticle to another 
(which can be not inevitable adjacent). Mathematically 
it can be proved in the following way. The processes of 
interaction in c (v) band can be considered to be indepen- 
dent on each other. It means, that transition probabil- 
ity from the (fc; )S (-state to (fe/ )S |-state in c-band and from 
(k'i J-state to (k'j s |-state in w-band, which is proportional 
to coefficient a2(k, k , s), can be expressed in the form of 
product of real parts of corresponding matrix elements, 
that is in the form 



a 2 (k,k',8) ~ Re(k^\V^\k^ s )Re(kl s \V^\k' 3}S ) 
Y,Re(KjV^\k ph )(k ph \k' n J, 



(11) 



where — Vo(„)e (e is unit operator) is the first term in 
Taylor expansion of pairwise interaction of quasiparticles, 
for instance, with wave vectors k' r , k' n and spin projec- 
tion s in w-band, that is, in ground state, V^°' — Vu c \ue 
is the second term in Taylor expansion of pairwise inter- 
action in excited state (in c-band), that is, it is product 
of configuration coordinate u and coordinate derivative 
at u — of operator of pairwise interaction of quasipar- 
ticles with wave vectors ki, kj and spin projection s in 
c-band, k p h is phonon wave vector, and the summation is 
realized over all the phonon spectrum. At that, since the 
linear density of pairwise interaction is independent on fc, 
which is the consequence of translation invariance of the 
chain, Vo( v \, Vu c ) are constants. Therefore, the pairwise 
interaction is considered to be accompanying by process 
of phonon generation, when electronic quasiparticles are 
already in excited state, that is, in c-band (retardation 
effect of phonon subsystem is taken into account) . Then 
we have — Vo( c \ue, — V ^e. A number of 

variants are possible along with process of phonon gener- 
ation, corresponding to states of electronic quasiparticles 
in c-band above described. The result will mathemat- 
ically be quite similar, if to change the energetic place 
of excitation, that is, if to interchange the role of c and 
v bands for given process. There seem to be possible 
the realization of both the stages (that is phonon gener- 
ation and absorption) for electronic quasiparticles in sin- 
gle c or v band states and simultaneous realisation both 
the stages in both the bands. Mathematical description 



will be for all possible variants similar and for distinct- 
ness, it was considered only the first variant. For the 
simplicity, the processes were considered, in which the 
spin projection is keeping to be the same. It is evident 
also, that in ^-direction the impulse distribution is quasi- 
continuous (since the chain has the macroscopic length 
L = Na). Consequently, standard way ^ fe h — > ^ f k ; 
can be used. Further, phonon states can be described by 
wave functions (k p h\ = voexp(ik p hz) , where z £ [0, L], 
k P h ^ [~ J ~' sH' v ° * s corLS t an t- Therefore, it is obtained 
from (jllf the expression 

a 2 (k, k',s) = b\vo v \ 2 \voc\ 2 V ( c )uVo( v )\(j)oca\ 2 \<f>ovs\ 2 x 
N 

— — -i?e{exp[i(fc;TO/ — kjmj)a] exp ika}x 

2ir(qi - qj){q r - q n ) 

Re{exp[i(k' r m r — k' n m n )a] expik'a}, 

(12) 

where |</>o cs | 2 , |</>ous| 2 are squares of the modules of the 
wave functions \kj }S ) and \k'j s ) respectively, k — k p h(qi — 
qj), k' = k' ph (q r - q n ) qi,qj,q r ,q n € N with additional 
conditions (qi — qj)a < L, (q r — q n )a < L, b - is aspect 
ratio, which in principle can be determined by compari- 
son with experiment. Here the values (qi — qj), (q r — q n ) 
determine the steps in pairwise interaction with phonon 
participation and they are considered to be fixed. We 
will consider The processes, for which k = k', are consid- 
ered. Consequently, (q r — q n ) = (qi — qj) and the operator 

^7r,u,a2 (^) is 



^7r,u,a2 (^0 — 

1212 12 4a2U sin ka sin k ' a ^' c) " +(v) " {v) " {c) 



c k',s c k,s C k,s> 



(13) 



where 4a2(s) is 



b | V 0v | 2 | V c | 2 V (c ) V ( V ) | 4>q cs | 2 1 (f>0v 

4a 2 (s) 



(14) 



2n[(qi - qj)] 2 
The system of operators c k , 



~+(v) Jv) -(c) 

c k>s > C U> c k's ^re- 



sponds to noninteracting quasiparticles, and it is under- 
standable, that in the case of interacting quasiparticles 
their linear combination has to be used 









-(c) 









ak,s 

Pk,s 



Otk,3 



-fc,s 

a(c) 



(15) 



Then the Hamiltonian ^ T f ai (u), which corresponds to 
SSH one-electron treatment of electron-phonon coupling, 
can be represented in the form 



^7r.u,ai (^) 



12 12 A *[^,« a M Ja M ~ a k,sPk, s al, 

k s 

+ /3 k . s a k ,sa k ; s >a k > s - P Ks a k ^ s >a k > s + a k s a k ^ s ' a\ ' s 

+ a k , s /3 k , s a k ^ 'a k ' s - /3 Ks a k , s a k ^ 'a\' s - P k , s a k ^ s 'a k> ' s ], 

(16) 



,2 a + (")a( c ) 
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where = Aaiusmka. The diagonal part V.f t ai (u) of 
operator 7^ W)U)Q:i (u) is 



^,t,a 1 («) = Z)Z)2A fc o MJ 9 M (i 



•(c) _ 



(17) 



where hjfl is density of operator of quasiparticles' num- 
ber in c-band, is density of operator of quasiparti- 
cles' number in w-band. The part of pairwise interaction, 
which is linear in displacement coordinate u and leads to 
participation of the phonons, is given by the Hamiltonian 

H (u) = ^^^4a 2 usinfcasinfc'ax 

k k' s 

(n 2 n +{c) n [v) - B 2 n {v) n +{c) 
\ a k',s a k',s a k',s Pk',s a k' : s a k',s 

X in? n + (c) n (v) - B 2 n Hv) n (c) 
x { a k,s a k,s a k,s Pk,s a k ,s a k,s 

The diagonal part r H. d K u Q2 (it) of operator H WjUj a 2 {u) is 
n d m («) = ±a 2 u J2 E E «* P* A (C) 

k k' s 

x otk,sPk,s{n^ s - n k c \) sinfc' a sin ka 



nw 



(19) 



Let us remark, that the Hamiltonian %TT,u,a 2 {u) de- 
scribes the attraction between the electrons, it can lead 
to formation of Cooper pairs in a 7r-subsystem and to 
superconductivity. 

The Hamiltonian - H 7Ti t(u) in terms of operator variables 

-(c) „(t>) 



2 u^+( c )a( c ) 



H„, t (u) = E E 2 *° cos fca K a L - Pl,s){at 



s "fc,s 



fc s 



a fe,J - 2 «fe,^fe, s (a fe a fe +a fc a. )j 



(20) 



The diagonal part t (it) of operator "Hjr.to (it) is given 
by the relation 



KtoM = EE e *(«L - - 4:1)> (2i) 

where = 2io cos ka. 

The operator transformation for the cr-subsystem, 
analogous to ([To]) shows, that the Hamiltonian Ho(u) is 
invariant under given transformation, that is, it can be 
represented in the form, given by ([7]). 

To find the values of elements of matrix in relation II 51) . 
the Hamiltonian H(u) has been tested for conditional 
extremum in dependence on the variables a k , fi k with 
condition s + 01 s = 1. 



Two values for the energy of quasiparticles, indicating 
on the possibility of formation of the quasiparticles of 
two kinds both in c-band and i>-band have been obtained. 
They are the following 



E ( k c \u) 



E [ k V \u) 



and 



V4 + Q 2 ^l 



(22) 



E { :\n) 



Q 2 H 



The factor Q is determined by relation 
QA k sin ka 



a 2 
2c7T 



where n^ s is eigenvalue of density operator of quasipar- 
ticles' number in c-band, ri^\ is eigenvalue of density 
operator of quasiparticles' number in w-band. The quasi- 
particles with the energy, determined by (|2U)) at Q = 1 
are the same quasiparticles, that were obtained in known 
SSH-model. 

Sufficient conditions for the minimum of the functions 
E(a kiS k e) can be obtained by standard way, which was 
used in [9j. It consist in that, that the second differ- 
ential of the energy being to be the function of three 
variables a ktSl p k ,s and X ktS has to be positively defined 
quadratic form. From the condition of positiveness of 
three principal minors of quadratic form coefficients the 
three sufficient conditions for the energy minimum have 
been obtained. Their analysis has showed, that SSH-like 
solution is inapplicable for the description of standard 
processes, passing near equilibrium state by any parame- 
ters. The quasiparticles, described by SSH-like solution, 
can be created the only in strongly nonequilibrium state 
with inverse population of the levels in c- and w-bands. 
At the same time the solution, the energy of quasipar- 
ticles for which is determined by (|22|) can be realised 
both in near equilibrium and in strongly non-equilibrium 
states of the 7r-subsystem of t-PA, which is considered to 
be quantum Fermi liquid. 

The continuum limit for the ground state energy of 
the t-PA chain with SSH-like quasiparticles will coincide 
with known solution [15|, if to replace A k Q — > A k . The 

calculation of the ground state energy Eq (u) of the t- 
PA chain with quasiparticles' branch, which is stable near 
equilibrium by taking into account, that in ground state 



EE 



(<0 



K, 



l k,s 



(23) 



)] = Q, (24) 



l k,s 



= 0, n v k — 1, in the continuum limit gives 



K\u) = - 



2Na 



(QA fc ) 2 4 dk + 2NRu 2 (25) 
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Then, calculation of the integral results in the expression 



E l o\u) 



4Nt Q 



TV 



{F(|,l-z 2 ) + 



2 „ w W 

- F(-,l- z 2 )}} + 2NKu 2 , 



where z 2 = 2^211* , F(§, 1 - z 2 ) is the complete elliptic 



2 _ 1Q.cl\u 
*° ' 

integral of the first kind and 1 — z 2 ) is the complete 
elliptic integral of the second kind. Approximation of 
ground state energy at z <C 1 for the stable near equilib- 
rium solution gives 



7T 



28(Q ai ) 2 u 2 
irt 



(27) 



+ ...} + 2NKu 2 . 



It is seen from ([27]). that the energy of quasiparticles, 
described by given solution, has the form of Coleman- 
Weinberg potential with two minima at the values of 
dimerization coordinate uq and —uq like to the energy 
of quasiparticles, described by SSH-solution. It is un- 
derstandable, that subsequent considerations, including 
electrically neutral S=l/2 soliton and electrically charged 
spinless soliton formation, that is the appearance of the 
phenomenon of spin-charge separation, by FL description 
of ID systems will be coinciding in its mathematical form 
with those ones in Fermi gas SSH-model. 

Thus, the possibility to describe the physical properties 
of ID systems in the frames of ID quantum FL includ- 
ing the mechanism of appearance of the most prominent 
feature of ID systems - the phenomenon of spin-charge 
separation - is proved. 

So, the model proposed takes into consideration the 
electron-electron correlations in explicit form, which 
seems to be ground for its application to electronic sys- 
tem of quasi- ID NTs, where electron-electron correla- 
tions are known to be rather strong. In particular, it 
can be used for analysis of ESR and FMSWR spectra. It 
can be done by above indicated manner, that is by using 
of hypercomplex number theory analogously to theoreti- 
cal analysis of quantum optics effects in Q and analysis 
in Q of Raman spectra in quasi- ID NTs. It seems to be 
also substantial, that FL treatment of electron-phonon 
interaction extends the applicability limits of SSH-model 
of ID conjugated conductors allowing its use in the case 
of strong electron-phonon interaction and (or) strong 
electron-photon interaction. It is evident, that the mech- 
anism of the phenomenon of spin-charge separation in ID 
FL is topological soliton mechanism, which is quite differ- 
ent from Anderson spinon-holon mechanism. Anderson 
mechanism of spin-charge separation is realised in fact 
for the systems with the linear function of the electronic 
energy in dependence on quasi-wave vector k for rather 
large energy interval. Just 2D NTs and graphene possess 
by given peculiarity, although the corresponding energy 
interval is relatively narrow. 



The results obtained show, that the shape of 7r-solitons 
(or er-solitons) is given by the expression with the same 
mathematical form both in Fermi gas SSH-model and in 
its FL consideration. It is 



W»)l s 



1 



£tt(<t) 



_ sech2r (n-n )a , , 



- ?v ((T) t]cos — , (28) 

S7T(ct) ^ 



where n, no are variable and fixed numbers of Cif-unit in 
CiJ-chain, a is C — C interatomic spacing projection on 
chain direction, Vn{o) is 7r (cr)-soliton velocity, t is time, 
^(a) is 7t(<t) coherence length. It is seen, that 7r-solitons 
(<7-solitons) differ in fact by numerical value of coherence 
length in Fermi gas SSH-model and in its FL description. 
Really, the coherence length £„• and £o<r is determined by 
the relation [22j 



Hvf c hvF 

S07T = ~7 j ?0(T = ~T j 



(29) 



where Ao CT) arc er— and 7r-bandgap values at T — 
OK, vp is Fermi velocity. In SSH Fermi gas model vp 
is proportional to to, in SSH FL-model, here presented, 
Dp is proportional to sum to + t\. It allows to explain 
some discrepancy between theoretical value in Fermi gas 
model and experimental values for £ w and its dispersion, 
depending on production technology. Theoretical value 
in Fermi gas model for £„. in t-PA is 7a, and it is low 
boundary in the range 7a — 11a, obtained for £ n from 
experiments [l8l |. It means, that in the samples with ir- 
soliton coherence length, equaled to 11a t\ is equal to 
0.57io at the same t-PA band gap value (it is possible, 
since factor Q is independent on ti and can be close to 
1). It means, that constant component t\ in Taylor ex- 
pansion of electron-electron interaction potential is sub- 
stantial and that it depends on preparation technology. 

The observation of FMSWR phenomenon in the NTs, 
formed by (110) high energy nickel ion implantation in 
natural diamond single crystal (at ion beam dose in 
5 x 10 14 cm -2 ) unambiguously means, that the SSH tt- 
soliton lattice formation takes place, that is periodical 
non-harmonic spin-density wave. For the simplicity we 
can approximate the function of the shape of given den- 
sity wave (let us designate it with abbreviation ASLSDW 
(Approximation of Soliton Lattice Spin Density Wave) 
by trapeziform function, FigfTJ] Real function of the en- 
velope shape of soliton lattice spin density wave will be 
between given trapeziform function and harmonic sinu- 
soidal function for density wave HSDW (Harmonic Spin 
Density Wave). We will show, that given approxima- 
tion is rather good for qualitative and in the princi- 
ple for quantitative estimation of characteristics of ferro- 
magnetic order formation. Actually, from mathematical 
viewpoint the difference between these states is minimal 
since ASLSDW can be represented in the form of super- 
position of HSDWs, where amplitude of main harmonic 
is strongly exceeding the others. It follows from the ex- 
pansion of ASLSDW profile in Fourier series, which is 
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f{x')=A(2L + L 1 ) 



OO 

E 

m=0 



[2(2L + L 1 )(cos 5 ^ Tr7 



cos ■ 



7rm(Lo+L 



2L +L 



7rLimsinm7r 



2Li7r 2 



cos ■ 



L 



(30) 



/(x) = A(2L„ + £i) 



(2L + £1 



E 

m— 1 



COS 



7rmLo 



2L +Li 



COS 



7rm(L +Li) 
2L +Li 



2Li7r 2 m 2 



27ttox 
L 



(31) 




Figure 15: Approximation of soliton lattice spin (charge) den- 
sity wave profile along the chain direction 



The sense of the parameters Lo, La, L is seen from 
FigCSl The functions f(x), f(x') denote spin density dis- 
tributions S(x), S(x') in the case of ASLSDW and charge 
density distributions p(x), p(x') in the case of ASLCDW. 
Expression (|30[) corresponds to position of chain atoms 
along AT '-direction, FigHU Expression (|3"Tj) corresponds 
to position of chain atoms along AT-direction, FigHSl We 
see, that by asymmetric deviation of spin (or charge) den- 
sity relatively the chain direction, the constant compo- 
nent is appeared. It will lead in the case of real SLSDW 
to ferromagnetic ordering if density wave distribution is 
magnetic spin density wave distribution and to ferroelec- 
tric ordering by similar electric dipole density distribu- 
tion. One can also see from (|30|) and ([3"Tj) . that the ampli- 
tude of components with m > 1 is inversely proportional 
to to 2 , that is, it diminishes rather quickly when num- 
ber to increases. In the first order approximation, the 
relations (|30p and (|5T|) can be represented by 

f(x')=A(2L + L 1 )acos—-+o(^), (32) 

/(*) = A (2L + LO (I + ~a cos + o (-Lj , 

(33) 

where a is 



It is seen from (|28p , that envelope function of single soli- 
ton spin density distribution corresponds precisely to the 



case with x axis position, that is, it is strongly asym- 
metric relatively the chain axis. By the formation of the 
soliton lattice the periodic envelope function of soliton 
spin density distribution in soliton lattice will be deter- 
mined by overlapping of individual envelope functions 
and the asymmetry extent can be even more pronounced 
depending on overlapping extent. Let us also remark, 
that given conclusion is preserved, when simultaneously 
with Peierls metal-insulator transition the Mott-Hubbard 
metal-insulator transition also takes place, it leads to 
modification of oscillating factor in (|2"5|) . Instead of factor 
cos 2 Tip the factor (b\ cos 2 ^ — b 2 sin 2 ^-) is appeared, at 
that in real case of t-PA b\ is substantially exceeds 62, al- 
though it is nonzero [23| . Further, it was found also, that 
there is numerical coincidence of the characteristics of 
SSH topological 7r-solitons in t-PA and in quasi-lD CZS- 
NTS, formed both by (111) and (110) ion implantation. 
So, the values of g-tensor components in Cu-implanted 
sample are g\ = 2.00255 (it is minimal g-value and it is 
<7|i principal direction), g 2 = 33 = g± = 2.00273, the ac- 
curacy of relative g- value measurements is ±0.00002 
g-value of paramagnetic 7r-solitons in t-PA, equaled to 
2.00263 (l9j . gets in the middle of given rather narrow 
interval of <?-value variation of PC in ion produced NTs. 
Although anisotropy of paramagnetic 7r-solitons in t-PA, 
mapping the distribution of 7r-electron density along the 
whole individual t-PA chain, is not resolved by ESR mea- 
surements directly, there are indirect evidences on axial 
symmetry of ESR absorption spectra in <-PA too, [20| . 
[21] • Consequently, the value 2.00263 is mean value and 
it coincides with accuracy 0.00002 with average value of 
aforecited principal g-tensor values of PC in quasi- ID 
CZSNTs. Taking into account given high precision co- 
incidence of g-values in t-PA and quasi-lD CZSNTs we 
can conclude that SSH-7r-soliton density distribution in 
quasi-lD CZSNTs is also asymmetric along hypercom- 
plex x-axis. Then, the appearance of ferromagnetic or- 
dering in quasi-lD CZSNTs, formed by (111) and (110) 
ion implantation is explained by a natural way. 

Therefore we have found the physical origin of the 
mechanism of the formation of ferromagnetic ordering in 
carbon NTs. It determines the only by 7r-subsystem of 
full electronic system, which agrees well with experimen- 
tal results on FMSWR above summarised and explains 
qualitatively the sensitivity of FMSWR parameters, pre- 
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sented in Figures 1 to 4 to symmetry directions of sur- 
rounding diamond crystal lattice by (110) implantation. 

Experimental results obtained by ESR study of NTs, 
formed in diamond single crystals in the result of the 
(100) ion beam modification indicates, that for the cor- 
rect description of Nts' properties the a-electronic sub- 
system has to be taken into consideration. It follows 
immediately from the appearance of inherent symme- 
try directions in distinction to (110) ion beam modifica- 
tion, whre ferromagnetic ordering symmetry characteris- 
tics are governed by symmetry directions of surrounding 
diamond lattice. 

Let us consider the following Hamiltonian 

H(u) =h (u) +n x (u) +n 2 (u) +h s (u) +n 4 (u), (35) 

where u is configuration coordinate, which corresponds 
to translation of CH-groups along the symmetry axis z 
of the chain, it is suggested to be independent on site 
position and on subsystem kind. Operator "Ho (it) is 

^w=EEEE £ ""Wl s %, ( 36 ) 

l rn q s 

in which subscripts m,q — {tt, <r}, s is spin projection, k 
is wave vector, ci , cr are operators of creation and 
annihilation of the quasiparticle with spin projection s 
and wave vector k in mth (gth) subsystem, e mq (u) are 
the resonance interaction integrals (hopping interaction 
in tight-binding model approximation) of quasiparticles 
in 7r-subsystem of t-PA electronic system, in cr-subsystem 
of t-PA electronic system which are considered to be 
Fermi liquids, and between tt- and a-subsystems. 
i-L\ is 

= EE Ul U> m > U ) d t mS (tf^(t) d trnsafjrns(i), 
m j 

. (37) 

where j = 1, N is the site position U\ (j, m, u) is intrasub- 
system Coulomb coupling parameter, which is dependent 
in general case on j, m, u. The operator 7-^2 («) is 

n 2 (u) = EE U ^ m > & M ) E £ tm S Cjms E ^qs^qs 
m>q j s s 

(38) 

where U%(j, m, q, u) is intersubsystem Coulomb cou- 
pling parameter, which is dependent in general case on 
j, m, q, u. The operator H^iu) is 

fisiu) = YYYYJi{i,m,q,u)c+ ms c+ qs ,c jms ,c 3qsi 

rn>q j s s 1 

(39) 

where Ji(j,m,q,u) is the inter-subsystem Hund's rule 
coupling, which is dependent in general case on j, m, q, u. 
The operator T-L^{u) is 

^4 = E E J 2(i' m '?' U )^ ms (t) 5 jm S (;)^a) 2 J9«(t) 
m^q j 

(40) 



where J2Q, m, is pair hopping parameter, which is 
dependent in general case on j, m, q, u. 

Like to foregoing consideration the Hamiltonians 
T-L\(u) and %2{u) can be expanded in Taylor series about 
the dimerizcd state. So, we have 

%(«) = EE(< , + 

m j (41) 
(-l) J 2a"V)c+ ms ^ t) Cj ms (^)C+ msW) Cj ms (^), 

where {a™}, m = {7r,er} are constants of electron- 
phonon interactions, accompanying the processes of in- 
trasubsystem Coulomb inter ations. 

M«) = E E( f/ 2 0) + 

(42) 

s s 

where {a™ 9 }, m, q — {tt, a} are constants of electron- 
phonon interactions, accompanying the processes of in- 
tersubsystem Coulomb interactions. It is clear, that the 
second terms in (|4"Tj) and in (|4"2"j) describe the attraction 
between strongly correlated electrons, it can explain the 
nature of the pairing mechanism in high temperature su- 
perconductors. 

The Hamiltonian T-L(u) can be considered to be basic 
Hamiltonian for description of the properties of carbon 
NTs, produced by (100) ion beam modification of dia- 
mond single crystals, in particular for analysis of ESR 
data above described. The generalization of the Hamil- 
tonian i-L(u) can be done by the way, proposed in Q on 
the basis of hypercomplex number theory, at that it has 
to be taken into account, that, strongly speaking, NTs, 
produced by (100) implantation cannot be described in 
the framework of ID theory, since C4 symmetry indicates 
on inequivalence of the chains along the NT axis. Let 
us remark, that the Hamiltonian (l35j) is similar to two- 
orbital-Hamiltonian, proposed in [27] for spectral analy- 
sis of the iron-based superconductors (SC). It will be al- 
most coinciding in the case when {a™} = 0, m = {tt, ct}, 
{a™ 9 } = 0, m,q = {7r,er}, Jx{j,m,q,u), J 2 (j,m,q,u) 
are independent on j,m,q,u. The difference in given 
case consists in inequivalence of a and tt- subsystems, in 
comparison with equivalence of Fe orbitals d xz and d yz , 
considered in [27^ . However even in given more simple 
case the task was solved the only by numerical methods. 
The main result is represented in Figure 8 in [271 ] . 

The magnetic excitation spectrum carries information 
on the nature of magnetism and the characteristics of 
superconductivity. It has been discussed in the litera- 
ture that an observation of a sharp quasiparticle-like res- 
onance peak in the spin fluctuation spectrum with the 
onset of superconductivity may strongly indicate a sign 
change in the gap structure caused by the superconduct- 
ing coherence factors. It has been established, that in 
iron pnictides a strong spin resonance occurs in the s- 
wave SC state. There is indication, that the spin reso- 
nance phenomenon is compatible with the s ± -wave SC 
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gap. The comparison of the ESR spectra with theoreti- 
cally calculated spectral function, presented in Figure 8 
in [23| allows to suggest, that the spontaneous transition 
in ESR response in the sample studied indicates on tran- 
sition to SC state, which is coexisting with antiferromag- 
netic (AFM) ordering. Therefore two lines - L-line and 
i?h-line - can be assigned with AFM-resonace observed 
in SC state. Inequivalence of the main characteristics of 
given lines can be attributed to strong inequivalence of 
two subsystems in NTs in comparison with theoretically 
calculated in [27|. It seems to be essential the conclu- 
sion in given work on the appearance of peak-dip-hump 
feature in a spectral function when the system becomes 
supcrcononducting. We have observed two dips. It seems 
to be consequence of different coupling of the resonance 
mode to fermions in it and cr-subsystems. Further the 
authors of [27j write " detailed study of the magnetic and 
the electronic spectrum shows that the dispersion of the 
magnetic resonance mode in the nearly isotropic su- 
perconducting state exhibits anisotropic propagating be- 
havior in an upward pattern" . Given conclusion is also in 
agreement with experiment. The observation of superlin- 
ear dependence in absorption kinetics is strong evidence 
of the mobility of spin carriers. The switching from the 
saturating behavior to nonsaturating behavior with su- 
perlinear absorption kinetics of -Rb-line can be attributed 
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